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Generating squeezing in an atom laser through self-interaction 
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We describe a scheme for creating quadrature- and intensity-squeezed atom lasers that do not 
require squeezed light as an input. The beam becomes squeezed due to nonlinear interactions 
between the atoms in the beam in an analogue to optical Kerr squeezing. We develop an analytic 
model of the process which we compare to a detailed stochastic simulation of the system using phase 
space methods. Finally we show that significant squeezing can be obtained in an experimentally 
realistic system and suggest ways of increasing the tunability of the squeezing. 
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Introduction. — The creation of the optical laser and 
the development of quantum optics has allowed tests 
of many fundamental properties of quantum mechanics 
UMM- The ability to create quadrature squeezing is an 
important prerequisite for many of these tests as it allows 
the creation of continuous variable entanglement between 
the amplitude and phase of two spatially separated opti- 
cal beams [2, 4]. With the advent of the atom laser, there 
is much interest in creating a quadrature-squeezed atomic 
beam as it allows us to revisit many of these tests using 
massive particles rather than photons. For example, it 
has been shown that continuous variable entanglement 
between the amplitude and phase of spatially separated 
atomic beams for use in Einstein-Podolsky-Rosen (EPR) 
tests [f| can be generated by dissociation of a molecular 
Bose-Einstein condensate (BEC) Q, or by outcoupling 
from a BEC using a Raman transition with squeezed light 
Q ■ Another example is interferometry — the use of mas- 
sive particles over photons already offers the promise of 
vastly improved sensitivity 0], and quadrature squeez- 
ing offers the possibility of going beyond this to beat the 
standard quantum shot-noise limit [9(. 

The standard scheme to create a squeezed atom laser 
is to use a squeezed optical field to couple atoms out of 
the BEC and into the atom laser beam, attempting to 
transfer the quantum state of the light onto the atoms 
Tot 11, Such a scheme is challenging, as it requires 
squeezed light at the relevant transition frequencies of 
the atomic species making up the BEC. Obtaining use- 
ful amounts of squeezing at these frequencies is a hard 
problem, although recently there has been some success 
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In this Letter we describe a scheme to generate a 
quadrature-squeezed atom laser without a squeezed opti- 
cal field, thus removing a significant source of complexity, 
and model the effect for an experimentally realistic sys- 
tem. Our scheme utilizes the nonlinear interaction caused 
by atom-atom scattering to create a Kerr squeezing effect 
[171 ]. The rate at which the beam squeezes is dependent 
purely on the local density of the beam. As the atoms 
fall under gravity, the beam become more dilute, leading 
to a continuous reduction in the strength of the Kerr ef- 



fect along the length of the beam. This ensures the Kerr 
effect acts only for a finite time, preventing the squeezing 
from becoming degraded in the long term limit. 

The possibility of nonlinearities resulting in quadrature 
squeezing has previously been considered by Jing et al. 
using a zero-dimensional, single- mode analysis [l8[ , who 
found very little squeezing within the range of validity 
of their linearized model. It has not been considered 
for a realistic atom laser, taking into account multimode 
effects, non-Markovian behavior (i.e. back coupling of 
the beam into the BEC) and mode matching. 

The structure of this Letter is as follows: We develop 
a single-mode model to obtain an analytic expression for 
the squeezing. We then create a realistic model of a Ra- 
man atom laser and simulate it using a stochastic phase 
space approach. These simulations are compared to a 
spatially integrated version of the analytic solution to 
determine its predictive power. Finally we consider the 
implications for current experimental atom lasers. 

Analytic model. — We first construct a single-mode 
model of Kerr squeezing that admits an analytic solu- 
tion. We use the Kerr Hamiltonian 
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where x is the strength of the nonlinearity and a describes 
a bosonic field. In the Fock basis, the evolution of a sys- 
tem governed by dJ) is described by \tp(t)) = J2 n c n (t)\n), 
with c n (t) = c„(0) exp[— i(nu + xn(n — l)/2h)t\. 



In order to examine squeezing in this system we define 
the standard quadrature operator X^ = e^a + e^^aJ , 
where <j) is the phase angle at which the measurement 
is carried out. The variances of the for all <f> are 
unity for a coherent state, and consequently a state is 
squeezed if the variance is less than one for a particular 
value of 4>. If we assume that our system is initially in 
a coherent state \tp(0)} = then defining var(A^) = 
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effective Hamiltonian for the system is 
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At any given time there is an optimum choice of <fi that 
gives the best squeezing. Plots of the minimum value of 
var(A^) over time for a variety of nonlinear interaction 
strengths are shown in Fig. [1] where we have chosen <f) to 
give the lowest possible variance of X*. The time taken 
to reach best squeezing scales inversely with \ while the 
minimum variance is given by ~ oT 2 ^ for a > 3. In 
this model arbitrarily good squeezing can be obtained 
provided the number of particles in the system can be 
arbitrarily large. 
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FIG. 1: Single-mode Kerr squeezing as a function of time 
for a nonlinearity \ an d particle number N = a 2 , x = 0.1ft, 
a = V1000 (solid); x = 0.04ft, a = yTOOO (dashed); x = 0.1ft, 
a = V500 (dash-dotted); % = 0.04ft, a = VE00 (dotted). 



Stochastic simulation of an atom laser. — We now de- 
velop a realistic, multimode and spatially extended de- 
scription of an atom laser. A mean field analysis using the 
Gross Pitaevskii equation (GPE) will not be adequate, as 
the GPE is incapable of examining the quantum statis- 
tics of the system, and thus cannot say anything about 
squeezing. We therefore model the system using stochas- 
tic phase space methods This involves finding the 
master equation for the system and then converting to a 
specific representation of the probability distribution to 
obtain a Fokker-Planck equation (FPE). This equation 
can then be treated as a set of stochastic partial differen- 
tial equations (PDEs) which can be solved numerically. 

Our model is based on a Raman atom laser [2(J [2l|, 
[22| . After adiabatically eliminating the excited state, the 
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where *i(r) and *2(r) describe the trapped and un- 
trapped matter fields respectively, Ay the single-photon 
detunings, 8 the two-photon detuning, Q13 and ^23 the 
Rabi frequencies of the two optical fields, ko = k2 — ki 
the momentum kick imparted to the outcoupled atoms 
(taken to be downward in our simulations), ui the har- 
monic trap frequency and f/y = Airh 2 aij /m, where ay- 
is the s-wave scattering length between atoms in states 
|*i) and |*j). As the matter fields are position depen- 
dent, Eq. describes the full multimode nature of the 
problem including non-Markovian effects. 

We will work in the Wigner representation, but ig- 
nore third and higher order derivatives in the FPE as 
these terms do not have a simple mapping to stochas- 
tic PDEs, and can be assumed to be negligible when the 
field has a high occupation number [23|. This truncated 
Wigner approximation (TWA) will eventually fail, but 
over the timescale of our simulations the TWA was in- 
distinguishable from the exact analytic solutions we had 
in the single mode case. In addition, quantities acces- 
sible by multimode GPE simulations also agreed with 
the TWA solutions over these timescales. The stochastic 
PDEs describing the system in the TWA are 
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where f2 = f^ 3 r2i2/Ai3 is the two-photon Raman Rabi 
frequency, ip\ and ip2 are the c-number stochastic vari- 
ables corresponding to the quantum operators for the 
trapped matter field and the atom laser beam respec- 
tively, and AV is the spacing of the grid on which the 
problem is to be numerically simulated. The terms in- 
versely proportional to AV compensate for the mean field 
of the vacuum, which is non-zero in the Wigner approach. 
As the FPE has no second order derivative term, there 
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are no explicit noise terms in the equations. Noise still 
enters the problem, however, as we must include the cor- 
rect noise in the initial conditions for Eqs. (0| and ((SJ). 
We chose this initial noise such that it corresponded to a 
coherent state. In all simulations parameters appropriate 
to a Rb Raman atom laser were chosen, i.e. a = 5.77nm, 
m = 1.44 x lCT 25 kg, k = 2 x H^m" 1 , Q = 50rads _1 , 
uj = 80rads _1 and a condensate with 5 x 10 5 atoms. 
As the simulations were carried out in one dimension, 
we assumed a cross-sectional area of 1.2 x KT n m 2 , and 
scaled Uij accordingly. The BEC nonlinearity Un was 
set to zero; this restriction will be discussed later. The 
stochastic equations (|4|) and ((5]) were solved numerically 
using the open source package XMDS [3]. 

Unlike the single-mode, zero-dimensional analytic 
model discussed earlier, the beam of an atom laser is an 
extended object, so we cannot talk of a particle number 
in the beam given by N = \a\ 2 as we could in the analytic 
case. In a beam the relevant quantity is the local density 
p. To quantify the squeezing in the multimode case, we 
define the amplitude and phase quadrature operators by 
picking a particular spatial mode of the quantum field. 
We define the amplitude and phase quadrature operators 
as X* = b + St, where b = f** L*(z)^ 2 (z) dz. Here, L{z) 
represents the spatial mode in which we are interested, 
and J x 3 \L(z)\ 2 dz — 1. Physically, the form of L(z) would 
be determined by the form of the local oscillator used in a 
homodyne measurement. The details of how the squeez- 
ing would be measured in an experiment will be consid- 
ered in the Discussion section. Maximum squeezing will 
be observed when L(z) best matches the spatial mode of 
the field in which the squeezing occurs. We chose L(z) 
as a plane wave with wavelength and frequency that best 
matched the atomic beam over the region Z\ — > Zq, i.e. 
L(z) = e *(kLZ-u L t+4>)^ with kL = ko _ U 22 pm/k h 2 , and 

u>L = hkj j /2m + U22p/fi. In Fig. [2] we plot the results of a 
stochastic simulation showing the variance of (lower solid 
trace) X^ and X^ +n / 2 (upper solid trace). <f> was chosen 
to minimise the variance of X^ (the squeezing quadra- 
ture), and hence maximise the variance oiX^I 2 (the 
anti-squeezing quadrature) . Initially the quadrature vari- 
ance is unity, as only vacuum is present. As the beam 
traverses the region, the overlap between the atom laser 
beam and L(z) becomes high, and squeezing is measured, 
reaching steady state shortly after the beam front has 
completely passed through the region. As the system 
reaches steady state mode matching is achieved, result- 
ing in the reduced variance signifying squeezing in the 
mode L(z). 

We now consider to what extent our simple ana- 
lytic model correctly predicts the squeezing. To com- 
pare the analytic model with the multimode simula- 
tions we choose \ = U22 J \L(z)\ 4 dz, a = ^/N where 
N = Q \L(z)\ 2 (\^ 2 (z)\ 2 - 1/2AV) dz, and then aver- 
age the var(A^) predicted by the analytic model over a 



period of time which corresponds to the time required for 
atoms to pass through the region [z\, z 2 \. The results are 
also shown in Fig. [2J which compares the best squeezing 
and antisqueezing predicted by the single-mode analytic 
model to the results of the stochastic simulations. The 
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FIG. 2: Variance of the quadratures (lower trace) and 
X^+V2 ( upper trace ) 

as a function of time for the atom laser 
beam in a region 20 pm long well below the condensate. The 
beam reaches the region 9 ms after outcoupling begins. The 
solid traces represent the results from the multimode stochas- 
tic simulation, and dashed traces represent our single-mode 
analytic model. After 12 ms steady state is reached in the 
region of interest. 



antisqueezing is very well predicted, but the single-mode 
model predicts squeezing almost two times better than is 
actually seen. This discrepancy is largely due to the diffi- 
culty of mode matching the multimode beam to the local 
oscillator. Any mode matching discrepancy will have a 
larger relative effect on the measured squeezing than the 
antisqueezing. 

Discussion. — Although our numerical simulations 
were carried out in ID, our scheme still functions for 
a real, 3D system. The main difference between ID 
and 3D is that 3D allows the beam to have a trans- 
verse mode structure. Theoretically this structure is ir- 
relevant as the squeezing rate is purely dependent on 
the particle density per mode, independent of dimen- 
sion, and the atom laser tends to single mode operation 
in the long time limit. It is still desirable to have as 
little transverse structure as possible, however, as this 
reduces mode matching problems. We now consider a 
realistic, 3D example. A Rb Raman atom laser, such as 
the one described in [22j], with mean trapping frequency 
Cu = 27r(60 x 600 x 6OO) 1 / 3 , has an atomic density just be- 
low the condensate of po ~ 3 x 10 18 m~ 3 if a two-photon 
Rabi frequency fi = 5001'ads^ 1 is chosen. A Raman 



atom laser is minimally divergent 22j, so density scales 
only due to acceleration by gravity. After falling a dis- 
tance z the beam density is p — po / (l+m,y/2gz / hko) . As- 
suming fco = 3.2 x 10 7 m _1 and that the mode match re- 
gion is a section of the beam 2bpm in vertical extent 1 cm 
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below the condensate, there are ~ 1100 atoms in this re- 
gion. Using these numbers our integrated analytic model 
predicts var(X^) = 0.143, and var(X* +7r / 2 ) = 7.11, 
where we use a time-dependent \ 1° model the density 
decrease as the atoms fall. Using Bragg diffraction as a 
beam splitter, squeezing of this level leads to entangle- 
ment under the Reid-Drummond criterion [B[ . While the 
measured amount of squeezing will not reach this due to 
mode matching difficulties, it indicates our scheme is cer- 
tainly feasible. If the nonlinearity or the density of the 
beam could be increased, the squeezing would further in- 
crease. Possible mechanisms to accomplish this might be 
the use of Feshbach resonances to increase the nonlinear- 
ity , or the use of far-detuned light fields to focus the 
atom laser beam and increase the atomic density. 

The flexibility of this scheme is clear: As the best 
squeezing depends only on the local density of the beam 
and how long atoms have been in the beam when they 
are measured, and since it is possible to tune the outcou- 
pling strength, momentum kick and place of measure- 
ment, there is a large parameter regime over which good 
squeezing can be obtained. 

Our scheme relies on the output beam starting in a co- 
herent state. As the outcoupling process functions as a 
beam splitter some of the quantum statistics of the BEC 
will be copied onto the beam. Assuming the BEC itself 
begins in a coherent state, it will also exhibit Kerr squeez- 
ing. However, due to the BEC's much higher density, it 
will reach peak squeezing more quickly than the beam, af- 
ter which the squeezing will degrade as the nonlincaritics 
cause its phase to become uncertain. The long-time limit 
of such a process is var(A^f ec ) = 2N for any quadrature, 
where ./V is the BEC particle number. As the outcoupling 
is weak, the beam will only weakly reflect the quantum 
statistics of the BEC, but due to the high variances this 
could still degrade the squeezing of the beam. 

There are at least two ways to obviate this problem. 
The first is to reduce the nonlinearities in the condensate 
using a technique such as Feshbach resonances. The non- 
linearities need to be suppressed such that the minimum 
shown in Fig.[T]occurs at time comparable to the length of 
the experiment, meaning the suppression factor can eas- 
ily be extracted from Eq. ([2|). For example, in the case 
of the Rb laser described earlier, the BEC nonlinearity 
need be reduced by a factor of approximately eight hun- 
dred. The second approach is to ensure the condensate 
remains near a coherent state due to continuous measure- 
ment and quantum back action. For example, one could 
use the scheme described in Ref. [26] , where a weak light 
beam continuously measures the condensate density. 

A homodyne measurement is required in order to ob- 
serve quadrature squeezing [2JJ. We note that this may 
prove challenging as obtaining a strong local oscillator 
which docs not itself undergo Kerr squeezing may be 
difficult. A similar problem exists in detecting optical 
Kerr squeezing generated in nonlinear fibers, and can be 



solved by using an asymmetric Sagnac interferometer to 
slightly rotate the quadrature axis of best squeezing so 
that it lies along the amplitude quadrature [28| . Now no 
local oscillator is necessary as the squeezing appears in 
the beam intensity, and can be measured simply by per- 
forming measurements of the flux. The analogy for an 
atomic beam involves mixing the squeezed beam with a 
weak reference beam sourced from the same BEC, with 
the phase chosen to rotate the axis of best squeezing onto 
the amplitude quadrature. This could be achieved by 
outcoupling atom laser beams in two separate internal 
states (for example, the F = 1, tuf = 0, and F = 2, 
mp = ground states of 87 Rb) which spatially overlap, 
and then recombine them using a microwave transition. 
The relative phase of the two field can be adjusted by 
altering the phase of the outcoupling field. As the fields 
spatially overlap, the mode-matching will be automatic. 
We have performed an analytic, single-mode analysis of 
the interference of two Kerr-squeezed atom laser beams 
with different intensities derived from the same conden- 
sate. In the case where the relative intensity between 
the two beams was 0.5, and using the same parameters 
used to generate Fig. [2l we found an intensity noise of 
0.17 as compared to unity for a coherent state when the 
two beams were interfered with the appropriate phase. 
The intensity squeezing obtained by this method is ro- 
bust to changes in the intensity of the two beams with 
best squeezing obtained for beam intensity ratios between 
0.25 and 0.50. This intensity squeezing is of considerable 
interest as one of the prime applications for atom lasers 
is precision interferometric measurements which are ulti- 
mately limited by shot noise. 
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